Beyond Gaussian : A Comment 
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In this paper we treat a non-Gaussian integral and give a fundamental formula in 
terms of discriminant. We also present some related problems. 

This is a comment paper to arXiv:0903.2595 [math-ph] by Morozov and Shakirov. 
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1 Introduction 

Gaussian plays a fundamental role in Mathematics (including Statistics) and Physics and 
et al. We want to (or we should) overcome its high wall in this century, and so a try is 
introduced. 

In the paper pQ the following "formula" is listed : 

j J e-( ax3+bx2y+cxy2+dy3 )dxdy = —L= (1) 
where D is the discriminant given by 

D = — (27a 2 d 2 + 4ac 3 - 18abcd - b 2 c 2 + 4b 3 d) (2) 
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of the cubic equation 

ax 3 + bx 2 + cx + d = 0. (3) 

The equation (CQ) is of course non-Gaussian. However, if we consider it in the framework 
of real category then ([T]) is not correct because the left hand side diverges. In this paper we 
treat only real category, and so a, b, c, d, x, y are real numbers. 

Formally, by performing the change of variable x = tp, y = p for (TjQ) we have 

LHS of (ED) = / [ e- pS ( at3+bt2+ct+d )\p\dtdp 

r( at3+bt2+ct+d y\p\dp\dt 



\a\e a " do I — dt 

^/(at 3 + bt 2 + ct + d) 2 



by the change of variable a = \/ at 3 + bt 2 + ct + d p. 
Therefore we can conjecture that the formula may be 



' dx = -]== (4) 



/r {/{ax 3 + bx 2 + cx + d) 2 ' V^D 
under the change t — > x. Here C is a constant. 

In the paper we calculate the left hand side of (jlj) directly. 

2 Fundamental Formula 

Before stating the result let us make some preparations. The Gamma-function T(p) is 
defined by 

POO 

T{p) = / e- x x v - x dx (p > 0) (5) 



o 



and the Beta-function B(p, q) is 

B(p,q)= I x^il-xY^dx {p,q>0). (6) 
Jo 

Note that the Beta-function is rewritten as 



oo x p-i 



Bip.q) = I — ——dx. 

! (1 + x)p+i 



2 



See [2] in more detail. Now we are in a position to state the result. 



Integral Formula 

(I) For D < 



where 



(II) For D > 



where 



r {/(ax 3 + bx 2 + cx + d) 2 \f—D 



v 2'6 y 



dx = 



■-£= ( 7 ) 



r yj (ax 3 + bx 2 + cx + d) 2 
1 Is 



(III) C_ and C + are related to C + = y/3C- by the equation 

^(5,5) = r2B(\, I). (9) 

Our result shows that the integral depends on the sign of D, and so our question is as 
follows. 

Problem Can the result be derived from the method developed in [1] ? 

A comment is in order. If we treat the Gaussian case (: e -( ax2 + bx v+ c y 2 )^ then the integral 
is reduced to 

1 dx = —^= (10) 



l R ax 2 + bx + c y/— D 
if a > and D = b 2 - Aac < 0. Noting 

i ri v 2'2' 



(JTU|) should be read as 

f I 

J R ax 2 + I 



-dx= 2 -^$S 
bx + c yf^D 
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3 Discriminant 



In this section we make some comments on the discriminant ([2]). See [3] in more detail ([3j 
is strongly recommended). 
For the equations 



f(x) = ax 3 + bx 2 + cx + d, f'(x) = 3ax 2 + 2bx + c 

the resultant R(f, /') of / and /' is given by 

a b c d 

a b c d 

R(f, f) = 3a 2b c 

3a 2b c 

3a 2b c 

It is easy to calculate (TT2|) and the result becomes 

-R(f, f) = 21a 2 d 2 + 4ac 3 - 18abcd - b 2 c 2 + 4b 3 d = -D 



(12) 



(13) 



On the other hand, if a, 0, 7 are three solutions of /(#) = in (jlip . then the following 
relations are well-known. 

' a + / ? + 7 = _i 

a(3 + 07 + £7 = £ (14) 
I «/87 = -I 



If we set 



the discriminant Z) is given by 



A = (a-/3)(a- 7 )09-7) 



D = a 4 A 2 . 

Let us calculate A 2 directly. For the Vandermonde matrix 

/ 1 1 1 \ 



(15) 



(16) 



V 



a (3 7 

4 



-A 



we obtain 



A 2 



(-\V\) 2 = \V\\V T \ = \VV T 



3 a + (3 + 7 a 2 + (3 2 + 7 2 

a + f3 + 7 a 2 + /3 2 + 7 2 a 3 + /5 3 + 7 3 
a 2 + /? 2 + 7 2 a 3 + /? 3 + 7 3 a 4 + /3 4 + 7 4 



a 



1 -1 

^ IT 



{ (6c - 9ad) 2 - 4(6 2 - 3ac) (c 2 - 3 



bd)} . 



(17) 



This result is very suggestive. In fact, from the cubic equation 



ax 3 + bx 2 + cx + d = 



we have three data 

A = b 2 - 3ac, B = be - 9ad, C = c 2 - 3bd 
, and so if we consider the quadratic equation 

AX 2 + BX + C = 
then the discriminant is just B 2 — 4AC. This is very interesting. 
Problem Make the meaning clear ! 

4 Concluding Remarks 

In the paper we calculated the non-Gaussian integral (j3J) in a direct manner. Details of 
calculation will be published in [3]. 

In this stage we can consider the general case. For the general equation 



f(x) = a x n + a\X n 1 + • • • + a n _ix + a, 



n 



(18) 



the (non-Gaussian) integral becomes 




(19) 
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The discriminant D of the equation f(x) = is given by the resultant R(f, /') of / and 
/' like 

I n(n — 1 'I n. T 77. — ~\ \ 

(20) 



I , nln— 1) 

a 



^ = (-l)^W,/')/«o 



where 



R(f, /') 



a ai 
a 



a\ ■ ■ ■ a n _i a n 



nao (n — l)ai • • • a n -i 

nao (n — l)ai • • • a n -i 



nao (n — l)ai 



see (dU) and ([12]). 



For example, if n = 4 and n = 5 then we have 



Q O *3 Q 0/1 /I O OOO OOO /l 

D n= 4 = 256a a 4 — 4a 1 a 3 — 27a a 3 — 27a 1 a 4 — 128a a 2 a 4 + a 1 a 2 a 3 + 16aoa 2 a 4 

Q Q < O Q 00 00 00 

— 4aoa 2 a 3 ~ 4a 1 a 2 a4 + 144a a2a 3 a4 — 6aoa 1 a 3 a4 + 144aoa 1 a2a 4 — 192a aia3a 4 
+18000102^3 + 18a^a 2 a3a 4 — 80aoaia 2 a3Ci4, 



and 



6 



A A Q "i Q Q O O Q O O 

D n=5 = 3125a a 5 — 2500aQa 1 a 4 aj: — 3750aQa 2 a 3 aj: + 2000aQa 2 a 4 a 5 + 2250aQa 3 a 4 a 5 
— 1600a a 3 a 4 a 5 + 256a a 4 + 2000a a 1 a 3 a 5 — 50a o a 1 a 4 a 5 + 2250a aia 2 a 5 

o o '3 O OA 

— 2050a aia 2 a 3 a 4 a 5 + 160a aia 2 a 4 a 5 — 900a aia 3 a 5 + 1020a aia 3 a 4 a 5 — 192a ai<2 3 a 4 
— 900aQa 2 a 4 a5 + 825aQa 2 a 3 ag + 560alala 3 a1a 5 — \2%a^a\a\ — 630ala 2 ala4 : a 5 

q Q n c 0/10 Q Q 

+ 144a a 2 a 3 a 4 + 108a a 3 a 5 — 27a a 3 a 4 — 1600a a^a 2 a5 + 160a Oi a 3 a 4 a 5 

o 000 000 00 0/1 

— 36a a 1 a 4 a5 + 1020a a 1 a 2 a 4 a 5 + 560aQa 1 a 2 a 3 a 5 — 746a a 1 a 2 a 3 a 4 a 5 + 144a a 1 a 2 a 4 

2 3 2 2 3 32 32 22 

+24aoa!a 3 a 4 a5 — 6aoa 1 a 3 a 4 — 630aoOia 2 a3a 5 + 2Aaoaia 2 a 4 a^ ) + 356aoflia 2 a 3 a 4 a5 

— 80aoaia 2 a 3 a 4 — 72aoaia 2 a 3 as + 18aoai£i2& 3 « 4 + 108aoa 2 a5 — 72aoa 2 a 3 a 4 as + 16aoa 2 a 

+16a a 2 a 3 a 5 — 4a a 2 a 3 a 4 + 256a^ag — 192afa 2 a 4 a5 — 128a^a 3 ag + 144a^a 3 a 4 a 5 

A A '3 00 Q00 "30 Q Q Q/1 

— 27a x a 4 + 144a 1 a 2 a 3 a 5 — 6a 1 a 2 a 4 a 5 — 80a 1 a 2 a 3 a 4 a 5 + 18a 1 a 2 a 3 a 4 + IQa^^a^ 
—Aa\a\a\ — 21a 1 x a\ai ) + 18a^a 2 a 3 a 4 as — 4a^a 2 a 4 — 4a^a 2 a 3 a 5 + a^a^a^a^. 

However, to write down the general case explicitly is not easy. 

Problem Calculate (fT9l for n = 4 (and n — 5) directly. 

The wall called Gaussian is very high and not easy to overcome, and therefore hard work 
will be needed. 
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